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We have studied the probability distribution of the perimeter and the area of the fcth largest erased-loop in 
loop-erased random walks in two-dimensions for fc = 1 to 3. For a random walk of A' steps, for large A', the 
average value of the fcth largest perimeter and area scales as A^^'^ and A' respectively. The behavior of the scaled 
distribution functions is determined for very large and very small arguments. We have used exact enumeration 
for A' < 20 to determine the probability that no loop of size greater than £ is erased. We show that correlations 
between loops have to be taken into account to describe the average size of the A;th largest erased-loops. We 
propose a one-dimensional Levy walk model which takes care of these correlations. The simulations of this 
simpler model compare very well with the simulations of the original problem. 
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I. INTRODUCTION 

The classical theory of statistics of extremals deals with the 
distribution of the largest of many independent identically dis- 
tributed random variables [|l]^]. After some rescaling, this 
distribution tends to one of the three universal distribution 
functions, the so-called Gumbel distributions, independent of 
the details of the starting distribution of random variables |g]. 
The independence of random variables is, however, not a good 
approximation in many physical problems. The statistics of 
extremes of many random variables is relevant in many differ- 
ent physical contexts. In many of these it is important to take 
account of correlations, for example, in the study of earth- 
quakes iQ], weather records |Q], slow relaxation in glassy sys- 
tems [H], and persistence in random walks [^. In some spe- 
cial cases extremal statistics of strongly correlated variables 
can be determined exactly [pj. In general, however, the study 
of extremal distributions of correlated and strongly correlated 
random variables poses a rather non-trivial problem even in 
the simplest cases. 

This paper deals with the extremal statistics of variables 
with long-range power law correlations in the loop-erased ran- 
dom walks (LERW's) in two dimensions. Our interest in the 
LERW problem comes from the fact that it provides one of 
the simplest examples of self-organized critical systems. In 
the LERW problem, the length of the walk is first increased 
by one at each step, and then decreases by a random amount 
due to possible loop erasures. The size of erased loops has 
a power-law tail [^. This is, thus, similar to the sandpile 
model where one grain is added at each time step but the num- 
ber of grains leaving the pile has a power-law tail. Clearly, 
there are correlations in the sizes of erased loops at differ- 
ent times. These correlations are more pronounced for larger 
loops. Erasure of a large loop leads to significant decrease in 
the length of the erased walk, and hence a significant decrease 
in the probability of erasure of another large loop within a 
short time. We propose that the expected ratios of sizes of A:th 
largest loop with the largest loop is a good variable to quan- 
tify these strong coiTelations, and propose a one-dimensional 
Levy walk model which is then tested by simulations. 



The LERW problem was introduced by Lawler [ |10| ] as a 
more tractable variant of the self-avoiding walk problem. This 
problem is related to many well-studied problems in statistical 
physics: the classical graph-theoretical problem of spanning 
trees, the ^-state Potts model in the limit q ^ [O], and the 
Laplacian self-avoiding walk problem ||ll|]. Connection to the 
spanning trees also relates this problem to the abelian sandpile 
model of self-organized criticality [[l3||. Recently simulation 
of LERW has been used as a computationally efficient way 
to determine the dynamical exponent of the abelian sandpile 
model in three-dimensions [p4|]. The upper critical dimension 
of LERWs is known to be 4 [|l5]|. In two dimensions, the 
fractal dimension of LERWs is known to be 5/4 p^ , [lq , p^ , 
and the exponent characterizing the probability distribution of 
the area of erased-loops is known to be superuniversal [nj]. 
Several other results on LERWs can be found in lP,p^-|20[] 
and a good review of earlier results on the LERW problem 
can be found in [Oil. 

In this paper, we show that the asymptotic behavior of 
the probability distribution function Prob(£(*'|A^) that the A:th 
largest erased-loop perimeter in the first A^ steps has value i^''' 
is described by a fc-dependent scaling function with argument 
ii^) /N~/^, We determine the behavior of the scaling function 
for the largest loop for very large and very small values of its 
argument. A similar behavior is found for the loops ranked by 
the enclosed area, rather than by their perimeter. The proba- 
bility that there is no erased loop of length greater than a fixed 
value r varies exponentially with A^ for large A^. Enumerating 
all walks satisfying this property (for a fixed r) is a general- 
ization of the self-avoiding walk problem. We have used ex- 
act enumeration techniques to determine the behavior of this 
probability for r = 0, 2, and 4 by enumerating all random 
walks with A^ < 20. We have proposed a simple Levy walk 
model which captures correlations in the LERW and agrees 
well with its extremal statistics as determined from large-scale 
Monte Carlo simulations. 

The plan of this paper is as follows: In Sec. pi after defin- 
ing the LERW model precisely, we recall relevant points from 
scaling theory for distribution of sizes of erased-loops. These 
are used to get the scaling form for the probability distribution 



of the perimeter and the area of largest erased-loop in a walk 
of A^ steps. In Sec. |l|, we outline our results about the con- 
nectivity constants /J2 and /J4 and determine their numerical 
value using series expansions. The simulation technique and 
results obtained thereof are described in Sec. |^. In Sec. M 
we describe the Levy walk model and compare the results of 
numerical simulations of this model with that of the LERW. 
Finally, some concluding remarks follow in Sec. ^|. 



II. SCALING THEORY OF LOOP-SIZE DISTRIBUTIONS 

A loop-erased random walk is defined recursively as fol- 
lows: For a one step random walk, the corresponding loop- 
erased random walk is the same as the random walk. To 
form the LERW l' corresponding to a given random walk 
of (A^ + 1 ) steps, we first form the LERW L corresponding to 
the first A^ steps of the random walk. Let us say this LERW 
L has n steps. We now add the (A^ + l)th step of the random 
walk to i . If no loop is formed, the resulting n + 1 stepped 
walk is l' . If this results in forming a loop of perimeter I, this 
loop is erased, and the resulting « + 1 — £ stepped walk is l' . 
A simple example is depicted in Fig. H 

Let L he a LERW of n steps obtained from a random walk 
of A^ steps. For a fixed N, n is a random variable. The critical 
exponent z of the LERW is defined by the relation that 



{n)r^N' 



z/2 



(1) 



for large A^, where the angular brackets denote ensemble aver- 
aging over all random walks of A^ steps. Since the root-mean- 
square end to end distance R for LERWs is the same as that 
for random walks, we have R ^ N^'^, and (n) ^ R^. Thus, z is 
the fractal dimension of the LERW. 

Let F{i\N) be the cumulative probability that a loop of 
perimeter greater than £ will be erased at the A^th step of the 
loop-erased walk. We define 

F{i) = lim F{£\N). (2) 

It was shown in Q that large N ':$> i '^ I, F{i\N) satisfies 
the scaling form 



F{e\N)^r^/\f{i/N-/^) 



(3) 



The scaling function f{x) tends to a non-zero constant as x 
tends to zero, and decreases to zero exponentially fast for 
jc 3> L Note that the exponents appearing in this scaling form 
depend only on the fractal dimension z- 

Let <&(£('' < £\N) be the cumulative probability that the 
perimeter £''' of the largest erased-loop in an A^-step walk will 
be less than or equal to L We shall study the behavior of this 
function for large A^. The probability that the largest erased- 
loop at the ki\\ step of an A^-step walk has perimeter less than 
or equal to (. is given hy \ — F{l\k). A simple approximate 
formula for (^{l'^ "^ < i\N) is obtained by neglecting correla- 
tions among sizes of erased-loops, and treating the generation 
of loops at different time steps as independent events. In the 



following, we will denote by Ouc the value of <I>(^('^ < £\N) 
in this uncorrected approximation. This gives 

N 

4>(^(^) < ^|A^) ~ 4>uc(£*^' < e\N) ^ Yl[l-F{£\k)]. (4) 

k=l 

Let X = £/N~/^, x(i) = ^^^/N"-!^ and y = k/N be new scal- 
ing variables. In terms of these new variables, substitution of 
F{£\k) from Eq. (g) in Eq. (|), gives 

1 



<i>uc(x(^'<x|A^)=n 



1 



N ^ \yzl^ 



(5) 



For fixed x and large, we can evaluate this expression by 
taking logs, expanding in powers of (l/N), and keeping only 
the lowest order terms in (1/A^). With this we get 

ln4>uc(^*^^ <x\N)^ -x-^l^f{x), (6) 

where fix) = /q / {x/y''^) dy. It is easy to see that f{x) has 
the same qualitative behavior as f{x). In terms of F{£\k), this 
equation can be written as 

,(^(') < ^|A^) ~ exp [-NF{e\N)] , 



1'u 



(7) 



where 



Fm) = -yF{i\k). 



N 



(8) 



k=0 



For small x, ln<I>uc(-^''' < x\N) should vary as —x ^''. For 
large jc, f{x) is small, and 1 — <I>uc(jc''' < x\N) should vary as 

Eq. (M) is a good approximation to Eq. (H) so long as the 
higher order terms in (1/A^) can be neglected. It is easily 
seen that the neglected term is of order A^F^(^|A^), and hence 
the approximation is valid so long as £ ^ N^-''^. It will be 
seen from simulation results (see Sec. |v|) that our assumption 
about correlations being small is not too bad and that Eq. (Q), 
and consequently also Eqs. (Q) and (0), are reasonable approx- 
imations to the largest loop size distribution for all £. The de- 
viation of the correct value from Eq. (^ is largest if £ is very 

small, say equal to 0, 2, 4, It is important to understand 

the behavior of <I>(^(^' < £\N) in this case. This we do in the 
next section. 

Let m be the expected number of loops of perimeter greater 
than or equal to £ generated from a random walk of A^ steps. If 
there are no correlations between different loops, for m <C A^, 
the number of such loops generated in particular realization is 
a random variable, distributed according to the Poisson distri- 
bution: The probability that exactly k such loops are generated 
is e^'"m^/k\. This implies that the probability that less than k 
loops of size greater than £ are generated can be expressed in 
terms of the probability that no loop of size greater than £ is 
generated. Simple algebra gives 



OucC^**^' < £,A^) = exp(-m) 



k-i 

7=0 



where 



= ^\og{<P,,(£<^'^<£\N)). 



(9) 



(10) 



ni. DETERMINATION OF CONNECTIVITY CONSTANTS 

Let Cr{N) be the number of A^-step random walks in which 
no loop of size greater than r is formed. The case r = corre- 
sponds to self-avoiding walks. As the total number of random 
walks of A^-steps is 4^ on square lattice, we have 



4>(£('' < r\N) = 



cm 



For large A^ it is expected that [E3 



(11) 



(12) 



For large A', Hr tends to a constant independent of A^, which 
may be called the "r-th connectivity constant" . From Jensen 
and Guttmann [pTJ the value of /jq is known very precisely 
and we have estimated /J2 and /J4 using series expansion and 
exact enumeration (details follow). 

Now consider the case k ~2. In this case, the walk cannot 
form loops, except that it can go back along the path it has 
taken. The connectivity constant /jq can be interpreted as the 
average number of forward directions available for the next 
step in an A^-step self-avoiding walk for large A^. As return 
along the direction of last step is now allowed, in the first ap- 
proximation, we should have 



A'2 « A'o + 1 



(13) 



We also have the trivial inequality /Jr < /Jr+2 for all r. As r 
tends to infinity, /Jr tends to 4. 

We determined the numbers Cr{N) for N <2Q and for all r 
by exact enumeration. The enumeration results for r = 2 and 
4 are tabulated in Table |. We analyzed this data by fitting it 
to the extrapolation form 



C,-{N)=KoiJ^{N)Ny-' 

Ki {-If 
N A^r+1/2 



1 



K2 



El 

N 



(14) 



where the critical exponent y is expected to be independent of 
r and takes the self-avoiding walk value of 43/32 in two di- 
mensions [ P4 I and Kj are constants which depend on r. This 
form is similar to that used by Conway and Guttmann [E3] 
for analyzing 51 -term series of self-avoiding walks. We have 
reduced the number of parameters in Eq. ( p^ because our se- 
ries is shorter Our estimates of /J2 and /J4, by fitting the form 
given by Eq. ( p^ term-by-term to the 20-term series tabulated 
in Table | are 3.7083(2) and 3.8818(4), respectively. These 
values are not very sensitive to variation in the fitting values 
of the parameters Kj. 

It is interesting to compare the numerical values of /jq, 
112, and /J4 with the estimates obtained using the uncorre- 
cted approximation. From Eqs. (Ill]) and ( [l2| ) we see that 
4>(€(i) < (.\N) varies as (a'^/4)" for large A^. Thus the ap- 
proximation Eq. (0) gives j^k/^ ~ 1 —F{k). Using the values 
of fj I; determined above, this would imply that F{Q), F{2), 
and F(4) have the values 0.3404, 0.0729, and 0.0295, respec- 
tively. The values of these quantities obtained from simula- 
tions are 0.3125, 0.0625, and 0.0257, respectively. We see 



that the approximation fares rather well in relating the prop- 
erties of the self-avoiding and loop-erased walks, which have 
quite different large-scale properties. 



IV. COMPUTER SIMULATION RESULTS 

We generated two-dimensional loop-erased random walks 
using the algorithm outlined in [|l4|]. For each walk we col- 
lected statistics about the perimeter and the area of the erased- 
loop at each step. The statistics were collected for A^-step 
walks with A^ = 2^ r = 15, • • •, 20. We averaged over 4.7 x 10^ 
different realizations of the random walk. We were able to 
simulate the entire ensemble in about 93 hours on a Pentium- 
Ill 700 MHz machine using about 2.6 Mb RAM. 



A. Largest loop perimeter 

During the simulations we collected statistics for F{£\N), 
the average number of loops of perimeter £ formed from a ran- 
dom walk of A^ steps. For each walk we also determined the 
perimeter and area of the five largest loops formed, this is used 
to obtain the measured cumulative distribution <I>o(^'*^' < ^1-^)' 
of size of loops of rank k, with fc = 1 to 5. The subscript "o" 
here refers to "observed". To reduce noise, nearby £ values 
were binned together We used 30 bins per decade of data. 

In Fig. we have shown the plot for Probo(^^*^' |A^) versus 
£ the observed probability distributions for k = 1, 2, and 3 
for A^ = 22". In Fig. | we have plotted <i>{£S'''^ < £\N) versus 
£/£* for various values of A^ as found in the simulations, and 
compared it to the theoretical curve given by Eq. (0) ignor- 
ing correlations between loops. An excellent collapse is seen 
among curves for all the values of A^ when plotted against the 
scaling variable x — £/£*. From these figures it is clearly seen 
that for X > 1 the prediction of the uncorrected theory is quite 
good and indeed asymptotically exact. However, considerable 
departure is seen for smaller values of x, forx ^ 1. 

We see that the prediction of the cumulative distribution 
function by the uncorrected theory is consistently higher 
compared to the observed distribution throughout the range 
of variation of the scaling variable x. This shows the expected 
anti-correlation between occurrences of large loops. 

For small values of the scaling parameter x, the observed 
cumulative distribution function seems to behave like 



4>o(x(') < jc) - flexp(-te"2/z) 



(15) 



with fl = 2.2 ± 0.3 and b = 0.39 ± 0.02. The fit is shown in 
Fig. ^ For large x, I - 4>o(x'^' < x) is very nearly NF{£\N) 
which varies as 



.(1) 



1 - <I>o(x^'^ < x) - flexp(-te^'-') 



2/z. 



(16) 



with the numerical value of the parameters obtained by curve- 
fitting being a = 0.32±0.03 and b = 1.7±0.1, same as that 
obtained by analysis of the all-loops data. This fit is shown in 
Fig.|. 



B. Largest loop area 

During simulations we collected statistics for the area of the 
erased-loops also. This was sampled exactly as the perimeter 
data in the previous subsection. 

In Figs, pi and M we have shown the plots for Probo(A(*^' |A^) 
versus A for A^ = 2^" and <1>(A('^) < A\N) versus A/N for var- 
ious values of A^, for fe = 1 to 3. The format of presentation 
is as in the previous subsection. An excellent collapse is seen 
among the curves for various values of A^ when plotted against 
the scaling variable y = A/N. 

The departure between the observed behavior and predic- 
tion of the uncorrelated theory is also similar to that seen for 
the perimeter data in the previous subsection. It is clearly seen 
from these figures that for y>0.l the prediction of the uncor- 
related theory is quite good and seems to be asymptotically 
exact for large y. For y < 0.1 considerable departure is seen 
between observed behavior and uncorrelated prediction. As 
in the perimeter data, there is a systematic over-prediction by 
the uncorrelated theory. 

For small values of the scaling parameter y, the observed 
cumulative distribution function <I>o(3''''' < y) seems to be- 
have like exp(— fl/y) with a = 0.049 ±0.002. For large y, 
1 —<i>oiy^^' < y) varies as e\p{—by) with /? = 14± 1. 



C. Variation of loop sizes with rank 

It is clearly seen in Fig. that the probability distribution 
of £'*^' becomes sharper as k increases. In fact, if k is of or- 
der A^ (say k = A^/1000), it is easy to see that the distribution 
tends to a delta function for large A^. A more careful argument 
shows that if A: ^ ^-vu+i)^ then the distribution would tend 
to a delta function. We note that i'^''' varies as {N/ky'^ and 
the average number of erased loops with this perimeter varies 
as N/{£^'^'y^^'^'. For the distribution to have sharp peak at 
fS \ this number should be much greater than fluctuations in 
the expected number of loops with perimeter greater than £'*^' . 
The latter varies as k^'^. Simple algebra then gives the re- 
quired result. 

A similar argument for the area distributions shows that the 
positions of the peak for the k-th rank varies roughly as N/k 
and their width varies as N/k^'^. Furthermore, when k^N'^'^ 
the width of the distribution becomes exponentially small in 
N. 



D. Affect of correlations on the probability distribution 
functions for fcth largest erased-loop size 

In Fig. |, we have plotted 0(^P) < £|^) ^^^ 4,(^(3) < ^|^) 
versus <!>(£(') < l\N) for A^ — 2^" from the observed distri- 
butions. This is compared with what would be expected on 
the basis of uncorrelated approximation. Similar plots using 
area (instead of perimeter) data show similar trends, and are 



omitted here. From this figure, it is clearly seen that the pre- 
dicted and the observed distributions are quite close. The ac- 
tual curve always lies above the value calculated by neglecting 
anti-correlations present. 



V. MODELING CORRELATIONS 

Consider the time series {«,} with / = 1, 2, . . . , generated in 
a LERW simulation, where «, is number of steps in the LERW 
at time step /. This process can be modeled by a stochastic 
motion of point on a one dimensional lattice. As rii is always 
positive, the motion occurs in the half space x > 0. In a single 
time step, this point can move one step to the right (if no loop 
erasure occurs in the corresponding random walk), or several 
spaces to the left. Now suppose that the random walk is not 
accessible to observation, and only the time series {«,} is ob- 
served. While the original LERW, treated as a stochastic pro- 
cess is a Markov process, the projected process is clearly not 
Markovian. However, it may be approximated as a Markov 
process. 



A. One-dimensional Levy walk model 

The transition probabilities for this Markov process are eas- 
ily defined. We think of «, as the position of a random walker 
at time / on a one dimensional lattice. The walk begins at 
f = with the walker positioned at x = 0. At each subsequent 
time step, the walker takes one step to the right and then draws 
a non-negative integer random number I with the probability 

Prob(^), ^ = 0, 1, 2, We will assume that for large £, 

Prob(^) decreases as f^"^ with T > 1. If ^ is less than or equal 
to the current position x of the walker, then the walker takes 
I steps to the left; otherwise it stays put. This completes one 
step. Clearly, we have 






(17) 



To ensure that there is no overall drift in the model, we also 
assume that 



££Prob(£) = l. 



(18) 



Note that the I here corresponds to the erased-loop size in 
LERWs. In general, one can expect to improve comparison 
with the original LERW model by making the probability of 
backward I steps when the walker is at n equal to the con- 
ditional probability in the LERW problem that the next step 
leads to erasure of a loop of length i when the current length 
of walk is n. This is expected to be of the form 



Prob(^|«) = Prob(£h)/cutoff(^/n), 



(19) 



where /cutoff is a cutoff function which is strictly zero if its 
argument is greater than 1 . We make the simple choice that 
/cutoff is 1 if the argument is less than 1 . 



For our simulations, we made a particular choice of 
Prob(^). We assumed that it is given by 



r 1 



Prob(^) = < 



1 



1 



£« (^+1)'^ 



for 1< ^ < oo 



(20) 



l-^Prob(/t), forf = 0. 

k=l 



Note that with this choice, the no-drift condition given by 
Eq. dl^ is automatically guaranteed for any choice of a. Fur- 
thermore, one can generate this distribution numerically by 
using only two calls to the random number generator We take 
a random number u with uniform distribution between [0, 1], 
define m = [m^''"J, and then put £ ~ m with probability 1/m 
and i = with probability 1 — 1/m. In our simulations, we 
used a ~ 0.6, which corresponds to the value T = 2.6 of the 
exponent of the two-dimensional LERWs. 

The Master equation for the above process describing the 
evolution of the probability P{x,t) of the walker being at po- 
sition X at time t is written as 



P{x,t + l) 



: 52Prob(£)P(x-l 

1=0 



(21) 



For large times f, the width of the probability distribution 
P{x,t) increases to infinity. It is easy to see that the width 
must increase as f '/(i+«). We note that if the particle it at x, 
its expected displacement in the next time-step is positive, as 
jumps with displacement greater than x to the left are disal- 
lowed. The contribution of such terms to Eq. ( |l7| ) varies as 
x^^'^. This equation may schematically be written in the form 






®p, 



(22) 



where 'D denotes diffusion operator which, presumably, in- 
volves fractional derivatives. The resulting equation for the 
scaling function is nonlocal, and its analytical solution seems 
difficult. Simple dimensional analysis shows that t scales as 
x^^'. Hence the width of this distribution should scale as 
f i/(t-i) Furthermore, for large f, P{x,t) tends to the scaling 
form 



P{x,t)'. 



1 



fl/(X-l)'' \tU{^-l) 



(23) 



B. Results from the Levy walk model 

We numerically integrated the Master equation Eq. ( pi| ) in 
X > half space using the probability distribution for erased- 
loop sizes given by Eq. (20) and computed P{x,t). The inte- 
gration for walks having up to A^ = 2^^ steps required about 
80 hours of CPU time on a Pentium II 350 MHz machine 
using about 7 Mb RAM. We also simulated the Levy walk 
process for time steps up to A^ = 2^" for obtaining the statis- 
tics on erased-loop sizes and the ktii largest erased-loop size. 



The quantities were sampled along the same lines as for the 
LERWs discussed in Sec. |^. To reduce noise in the statis- 
tics, we averaged over a large ensemble consisting of 2 x 10^ 
different runs. The simulation of the entire ensemble required 
about 141 hrs of CPU time on a Pentium II 350 MHz machine 
using about 1 .5 Mb RAM. 

Scaling plots for the computed probability of finding the 
Levy walker at location x at time step A^, P{x,N), are shown 
in Fig. ^ In this figures we have plotted N^'^P{x,N) versus 
x/N^''^, for z = 5/4. The figure clearly shows that the ob- 
served behavior agrees well with the conjectured scaling form 
given by Eq. @. 

We also analyzed the distribution of kth largest loop sizes 
in simulation of this Levy walk model, and compared them 
with the corresponding distributions for the two-dimensional 
LERW model. We found that the deviations from the predic- 
tions of the uncorrected theory are much smaller in the case 
of the Levy walk model than in the original LERW. The plots 
are very similar to the Figs. g[ ^ and || and are being omitted 
here. 

In Fig. |l^ we have compared the probability distributions 
for the ^th largest erased-loop sizes from the Levy walk model 
with those from LERW. The figure clearly shows that the 
probability distributions obtained from the Levy walk model 
match very well with those from the LERW. 

A better quantitative estimate can be obtained by comparing 
the ratio R^, defined as 



/?, = (#-))/(^(i)), 



(24) 



where () denotes expectation value. 

The value of R^ as found in the simulations of the LERW 
was found to be 0.605, 0.463, 0.386, and 0.335 fork = 2 to 
5, respectively. The corresponding values in the simulation 
of the Levy walk model were 0.614, 0.474, 0.397, and 0.347, 
respectively. The corresponding values from the uncorrelated 
approximation would be k^''^, i.e., 0.648,0.503,0.420, and 
0.365, respectively. It is clear that the Levy walk model gives 
a much better estimate of these ratios than the uncorrelated 
approximation. 



VI. CONCLUDING REMARKS 

Our analysis above shows that the probability distribution 
of the largest erased-loops in LERWs is fairly well described 
by the simple approximation ignoring correlations between 
the sizes of different loops. However, the average values of 
ratios of v^' are not well described in this approximation. A 
simple model which takes care of a large part of these corre- 
lations is the Levy walk model introduced in this paper. In 
this model, one keeps information about the length of the 
LERW, but throws out all information about its shape. We 
have seen that this model reproduces the extremal statistics of 
the LERWs quite well. 

Secondly, we have exactly enumerated Cr{N) the number 
of A^ step LERWs in which loops of size less than or equal 
to r are erased. Using these we have determined /j^ the rth 



connectivity constant. The determination of /jq for various 
lattices has been a long-standing problem in lattice statistics. 
Higher r-values present interesting geometrical questions, and 
may be helpful in understanding the crossover from random 
walk to self-avoiding walk. 
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TABLE L Number of A'-step loop-erased random walks Q (A') in 
which the largest loop of perimeter I less than or equal to 2 and 4 are 
erased for A' = 1, ■ • ■, 20. 



A' 


CiiN) 


Ca{N) 


1 


4 


4 


2 


16 


16 


3 


64 


64 


4 


248 


256 


5 


976 


1024 


6 


3736 


4072 


7 


14536 


16248 


8 


55280 


64352 


9 


213336 


256120 


10 


808016 


1011504 


11 


3099456 


4016496 


12 


11706568 


15828968 


13 


44696992 


62727520 


14 


168475176 


246805224 


15 


640913784 


976340664 


16 


2411998168 


3836482296 


17 


9148925856 


15153764480 


18 


34387933200 


59482843856 


19 


130125970320 


234640138528 


20 


488603502672 


920216177360 




FIG. 1. An illustrative example of the loop-erasure proce- 
dure and some aspects related to perimeter and enclosed area 
of erased-loops in loop-erased random walks: The random walk 
a-b-c-d-i-b-e-f-e-g-h-g-i-j-k-1 of 52 steps starts at a, and ends at 1. 
The erased-loops are shown by thin lines and the loop-erased walk 
a-b-l-j-k-1 having 12 steps is shown by thick lines with sites on it 
marked by solid circles. Note that at the points i and k, while the 
random walk path intersects itself, the LERW encounters no inter- 
section as the loop b-c-k-d-i-b has already been erased. 



^ 1.5 



Q 




8x10" 



FIG. 2. The observed probability distributions for perimeter of the k\h largest erased-loop, k= 1,2, and 3, for two-dimensional LERW for 

^ = 220. 




FIG. 3. The cumulative probability distribution for perimeter of the kth largest erased-loop, k = I, 2, and 3, for different values of A' for 
two-dimensional LERW. Solid lines give the prediction of the uncorrelated theory and dashed lines with symbols give the numerically observed 
distributions. For £/£* > 1 the curves match well with ^{fM' < £\N) approaching unity very fast. Note the excellent collapse of the lines of 
the same type for all values of A' and k and also the systematic deviation (over prediction) of the uncorrelated theory from the numerically 
observed distribution. 
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FIG. 4. Variation of tlie cumulative probability distribution for perimeter of the largest erased-loop for small II for different values of A' 
for two-dimensional LERW. Solid line gives curve-fit corresponding to Eq. dla) and dashed lines with symbols give the numerically observed 
distributions. 




FIG. 5. Variation of the cumulative probability distribution for perimeter of the largest erased-loop for large £ for different values of N for 
two-dimensional LERW. Solid line gives curve-fit corresponding to Eq. (Ilq) and dashed lines with symbols give the numerically observed 
distributions. 
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FIG. 6. The observed probability distributions for area of the kth largest erased-loop, k = 1, 2, and 3, for two-dimensional LERW for 

yv = 220. 
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FIG. 7. The cumulative probability distribution for area of the klh largest erased-loop, fc = 1, 2, and 3, for different values of N for 
two-dimensional LERW. Solid lines give the prediction of the uncorrelated theory and dashed lines with symbols give the numerically observed 
distributions. For A/N > 0. 1 the curves match well with <I>(A") < A\N) approaching unity very fast. Note the excellent collapse of the lines 
of the same type for all values of A' and k and also the systematic deviation (over prediction) of the uncorrelated theory from the numerically 
observed distribution. 
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FIG. 8. Variation of the cumulative probability distribution for perimeter of the fcth largest erased-loop, k = 2 and 3, with that of the largest 
erased-loop for two-dimensional LERW. Dashed lines give the prediction by uncorrelated theory and solid lines give the behavior of the 
observed data. Here the curves are shown only for N = 2 . Curves for other values of N = 2'', r = 17, 18, 19, collapse indistinguishably with 
these curves. 
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FIG. 9. Scaling plots from numerical integration of the Master equation Eq. (Ell) for probability of finding the Levy walker at position x at 
time step A' versus x/N^' , z = 5/4, for A' = 2 , and 2'^. Good scaling and consequently good collapse of curves is seen. 
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FIG. 10. Observed probability distributions for size (perimeter for LERW) of the klh largest erased-loop for two-dimensional LERW (solid 
lines) and the Levy walk model (dashed lines) for A' = 2^". The extremal distributions for the Levy walk model have been rescaled by 
multiplying (dividing) the abscissa (ordinate) by a factor of 1.04. This rescaling makes the mean points of the distributions obtained from the 
Levy walk model coincide with those of the LERW. 
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